Nonlinear shot noise in mesoscopic diffusive normal-superconducting systems 
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We study differential shot noise in mesoscopic diffusive normal-superconducting (NS) heterostruc- 
tures at finite voltages where nonlinear effects due to the superconducting proximity effect arise. A 
numerical scattering-matrix approach is adopted. Through an NS contact, we observe that the shot 
noise shows a reentrant dependence on voltage due to the superconducting proximity effect but the 
differential Fano factor stays approximately constant. Furthermore, we consider differential shot 
noise in the structures where an insulating barrier is formed between normal and superconducting 
regions and calculate the differential Fano factor as a function of barrier height. 

PACS numbers: 74.40.-|-k, 74.50.-|-r 
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Shot noise is fluctuation of the current that is due 
to the discrete nature of the charge carriers. It is the 
only source of noise at zero temperature. Current noise 
contains information on the physics of transport phe- 
nomenon not contained in the conductance. A classical 
value S = 2e|/| = Spoisson for the noise, observed, e.g., in 
a vacuum diode, is obtained in the tunneling limit, i.e., 
when the transmission probabilities for open scattering 
channels are small and there are no correlations among 
the charge carriers. However, two effects may reduce the 
noise below its Poissonian value: inelastic scattering (not 
considered in this paper) and reduced noise in channels 
with finite transmission amplitudes.u In a phase-coherent 
conductor, all the transfer coefficients are not necessarely 
small. Instead, e.g., in a diffusive phase-coherent metallic 
wire the distribution function of the transmission coeffi- 
cients has a bimodal forma such that almost closed and 
almost open channels are preferred. This results in a 
shot noise that is one-third of the Poissonian noise inde- 
pendently of the sample-specific properties SVgh as the 
number of channels or the degree of disorder .H'B 

During the last decade the noise properties of meso- 
scopic conductors have been under intense study (for a 
review, see Ref . ^) . There has also been increasing inter- 
est to comprehend the interplay of phase coherence and 
superconducting proximity effect in mesoscopic physics. 
Recently the doubling of the shot noise in normal- 
superconducting heterojunctions, predicted in the linear 
regime in Ref. |^, has been verified experimentally.B While 
the theory of classical shot noise and the quantum me- 
chanical results in the linear regime have been discussed 
before, little attention has been paid to mesoscopic noise 
at finite voltages when the presence of superconductivity 
induces nonlinear behavior in the transport coefficients. 
In Ref. current correlations in very small hybrid NS 
structures at finite voltages were discussed. Eliminating 
the effects due to the finite size of the structure and fully 
taking into account the effects arising in a diffusive phase 
coherent sample, however, requires larger structures to 
be studied. In Ref. ^ a counting-field approach to the 
Keldysh Green's-function method was adopted to calcu- 
late numerically the statistics of current in a normal wire 
connected to normal and superconducting reservoirs at 



finite voltages. In this paper we use a well-established 
scattering-matrix approach to calculate the differential 
shot noise at finite voltages in the presence of the prox- 
imity effect. 

We find that in the presence of superconductivity, the 
differential shot noise follows the reentrance peak ob- 
served in conductance such that the differential Fano fac- 
tor remains approximately constant. In the NS structure 
the resulting differential Fano factor is twice its normal 
value. The Fano factor can be roughly interpreted to 
be the ratio of the effective charge-carrying unit and the 
unit charge. Hence the doubling of the shot noise is a 
signature of Cooper-pair transport in the NS junction. 

In the second part of the paper, we consider differ- 
ential shot noise in structures where an insulating tun- 
neling barrier separates the normal and superconduct- 
ing regions. In the tunneling limit differential conduc- 
tance essentially probes the density of states in the su- 
perconducting side except that at zero voltage and low 
temperature the reflectionless-tunneling effect increases 
the differential conductance. Reflectionless tunneling 
arises because of the quantum coherence of electrons and 
Andreev-reflected holes that travel along the same paths 
in opposite directions scattering several timeS|-|fxopi. the 
barrier and the disorder potential of the metal.&LjO We 
calculate the differential shot noise as a function of volt- 
age and predict that the reflectionless-tunneling effect is 
present not only in the differential conductance, but also 
in the differential shot noise, such that the differential 
Fano factor stays constant. 

We consider the zero-frequency shot-noise power in a 
two-lead system that is the a; = limit of the Fourier- 
transformed current-current correlation function 



S = 



dt([/(t)-(/(t))][/(0)-(/(0))]), (1) 



where the current operator I{t) may be expressed 
through a scattering matrix s and () denotes the quan- 
tum mechanical expectation value. In the NS junction 
where one reservoir is normal and the other is supercon- 
ducting the two-terminal differential shot noise at zero 
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temperature T takes the for 
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Here s'lf i^ii) is the submatrix of the scattering matrix 
referring to the reflection as a hole (electron) of an elec- 
tron incident in lead 1, evaluated aX E — eV . AtT^O 
the differential conductance of the NS junction at voltage 
V is given by 



G — GqTt §2^82^^ + S21S21 ' + 2s1ls'l\ 
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(3) 



where the unit conductance of a single spin-degenerate 
channel is denoted by Gq = 2e^/h. In the absence of 
normal transmission for voltages below A/e Eq. re- 
duces to 



1 dS 

eGo dV 



(4) 



and the conductance is directly proportional to the An- 
dreev reflection probability, 



Gns = 2GoTr[sJjs'if ], 



(5) 



the factor of two indicating the fact that Andreev reflec- 
tion creates a Cooper pair in the superconductor. 

We have studied two kinds of NS heterostructures: 
a phase-coherent normal-metal wire connected from the 
one end to a normal reservoir and from the other to a long 
superconductor, and an NIS structure including a tunnel- 
ing barrier between the normal wire and the supercon- 
ductor. In our calculations we adopt a scattering-matrix 
approach and apply a numerical decimation method to 
truncate |tiie Green's function of the two-dimensional 
structure.ll3 The disordered normal-metal structure is 
modeled by a tight-binding Hamiltonian with the site en- 
ergies varying at random within range [— ^w, ^w] . Here 
we choose w — where 7 is the nearest-neighbor cou- 
pling parameter. The calculated values of the observ- 
ables are averaged over several impurity configurations. 
The parameters of the structure are chosen such that the 
transport through the normal metal is diffusive, i.e., the 
mean free path / is much smaller than the length L of 
the structure, which on the other hand is much smaller 
than the localization length Nl, where N is the number 
of quantum channels {I L <^ Nl). The length scales of 
the structures in units of the lattice constant are depicted 
in the insets of Figs. |l| and ^ illustrating the scattering 
geometries of the NS heterojunctions. 

The influence of the proximity effect on current noise 
is most conveniently seen in the differential Fano factor 
{dS/dV)/2eG. In order to compare the electrical trans- 
port through normal and normal-superconducting struc- 
tures we calculate {dS/dV)/2eG in these two cases. The 
calculated differential conductances for normal (Gn) and 
NS (Gns) structures are plotted as functions of voltage 
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FIG. 1: Differential conductance as a function of voltage. Cir- 
cles, NS structure depicted in the inset; squares, N structure. 
The 95% confidence interval for the relative error is ±1%. 
The results in Figs. |l| - ^ were ensemble averaged over 600 
realizations. 
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FIG. 2: Differential shot noise as a function of voltage. Cir- 
cles, NS structure; squares, N structure. The 95% confidence 
interval for the relative error is ±1%. 



in Fig. 0. As expected, in NS conductance we observe 
the well-known nonlinear reentrant behaviour due to the 
presence of the superconducting proximity effect, i.e, Gns 
exhibits a maximum at energies of the order of a few 
Thoulcss energies Et — HD/L^. At this energy scale 
Gn remains constant. At zero voltage, normal shot noise 
is known to have a value iSn = jtS'poisson- In the nor- 
mal case there are no nonlinearities at this energy scale. 
Thus increasing voltage does not change this result and 
the differential shot noise dSj^/dV for the normal struc- 
ture remains constant as shown in Fig. ^ In the NS case, 
however, the differential shot noise dS]<ss/dV exhibits a 
similar reentrant effect as Gns ■ This is in agreement with 
the counting-field approach of Ref. |[ Combining these 
two results, in the N structure, we observe that for the 
differential Fano factor the result {dS^/dV)/2eG =1/3 
holds also for finite voltages, i.e., differential shot noise 
has the value one third of the Poisson value. In the 
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FIG. 3: Differential Fano factor {dS/dV)/2eG as a function 
of voltage. The 95% confidence interval for the relative error 
is ±1%. Circles, NS structure; squares, N structure. The 
small variation in the NS case is outside the error bars, but 
relatively smaller than the corresponding variations in G and 
dS/dV. 




FIG. 4: Normalized difi'erential conductance G/Gb as a func- 
tion of voltage in the NIS structure depicted in the lower 
inset {Gk/Gb = 10). Circles rfcattering-matrix approach; 
solid line, quasiclassical theorya. Upper inset: difi'erential 
conductance at the larger voltage scale. The 95% confidence 
interval for the relative error is ±1%. The results in Figs. 4 - 
6 were ensemble averaged over 300 realizations. 



NS structure differential shot noise follows the reentrance 
peak observed in differential conductance, such that the 
differential Fano factor {dS-[<s /dV)/2eG Ri 2/3 remains 
approximately constant (Fig. Hence, the differential 
Fano factor is twice the normal value reflecting the fact 
that in an NS junction the current essentially results from 
the uncorrelated transfer of Cooper pairs. 

In Fig. 1^ we note that at zero voltage, Gns is 2% be- 
low Gn. This is due to a weak-localization effect result- 
ing from the quantum iaterference between time-reversed 
paths of the electronsllirllJ yielding a different contribu- 
tion to Gns than to Gn. 

If the N wire is only weakly connected to the su- 
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FIG. 5: Normalized differential shot noise {dS/dV)/eGB as a 
function of voltage. The 95% confidence interval for the rela- 
tive error is ±1%. Inset: Differential shot noise (squares) and 
differential Fano factor (circles) at the larger voltage scale. 




FIG. 6: Differential Fano factor (dS/dV) /2eG in an NIS struc- 
ture at 1/ = as a function of resistance of the tunneling 
barrier Rb/Rn (measured in units of normal structure -Rn). 
The 95% confidence interval for the relative error is ±1%. 



perconductor, the differential conductance G probes the 
density of states of the superconductor. This is de- 
picted in the upper inset of Fig. ^ where the conduc- 
tance of a NIS structure is plotted as a function of volt- 
age. However, at zero voltage, the coherent interplay be- 
tween the Andreev-reflected and disorder-scattered elec- 
trons results in the reflectionless-tunneling effect increas- 
ing Gnis and creating a peak around V = (Fig. Q). 
The reflectionless-tunneling effect may arise when there 
is a superconductor on the other side of the tunneling 
barrier and multiple scatterings by the tunneling barrier 
and by the-.di.sorder potential in the normal-metal area 
take place.Mm 

As a test of our numerical results, we compare the 
obtained differential conductance G in the case of re- 
flectionless tunneling to the expressions derived from the 
quasiclassical theory of nonequilibrium, inhomogeneous 
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superconductivitytSta in the diffusive limit. To calcu- 
late the current for the present setup, we need to solve 
for the transverse distribution function (i.e., the symmet- 
ric part of the electron distribution function around the 
chemical potential of the superconductor) whose gradi- 
ent determines the quasiparticle current. In the limit 
rf, = Rb/Rn > 1, < A, we obtain for G at T = 
(for details, see Ref. 16 and the references therein) 



G = 



GN{sa\2^ + sinh2\/t;) 



Ar'^^Jv{cos'^ y/v + sinh ^Jv) -f sin 2^ -\- sinh 2^ 

(6) 

where v = cV/Et and Gn is the normal-state conduc- 
tance of the diffusive wire. This is also plotted as a func- 
tion of V in Fig. M. and agrees well with the scattering- 
matrix approach .Ej 

We have also calculated the differential shot noise as 
a function of voltage in the NIS structure (Fig. ^) us- 
ing the scattering-matrix method. We observe that also 
the differential shot noise increases at zero voltage due 
to reflectionless tunneling. The inset in Fig. |^ illustrates 
the differential shot noise and the differential Fano factor 
at larger voltages. At voltages slightly above zero cur- 
rent and noise are suppressed since there are no single- 
particle states in the superconductor and the Andreev 
reflection probability is proportional to the square of the 
tunneling probability. At V = A/e the differential shot 
noise follows the conductance peak. The interplay be- 
tween the Cooper-pair and single-electron transport is 
most clearly seen in the differential Fano factor. At volt- 
ages below A/e the Andreev reflection is the dominant 
charge transfer mechanism. Since the transmission prob- 
ability through the barrier is small, at low voltages we ob- 



tain a Fano factor which is twice the Poissonian value. As 
voltage approaches A/e, the normal transmission proba- 
bility increases and the system effectively behaves more 
like a normal conductor. Thus, the differential Fano fac- 
tor quickly decreases to a value near unity. In the inset 
of Fig. 5 the Fano factor at F = is below the values 
obtained at somewhat higher voltages. This is due to the 
fact that the chosen value for r?, = 10 is not strictly in 
the tunneling limit. 

In order to illustrate the crossover from an ideal inter- 
face to the tunneling limit we have studied reflectionless 
tunneling by calculating the Fano factor as a function 
of the interface resistance Rb at zero voltage. The resis- 
tance Risr of the normal structure gives the characteristic 
scale for Rb, thus we plot the differential Fano factor as 
a function of Rb/Rn (Fig. ^. At large values of Rb, 
differential Fano factor approaches a limiting value two. 

In conclusion, we have studied differential shot noise 
in normal-superconducting mesoscopic structures in the 
nonlinear regime at finite voltages. The superconducting 
proximity effect manifests itself as a well-known nonlinear 
reentrance behaviour in the conductance at the voltages 
of the order of a few Et/b. We have shown that also 
the shot noise exhibits a similar reentrance effect which 
keeps the differential Fano factor approximately constant 
as a function of the voltage. In the second part of the 
paper, we have considered a nonideal NS interface with 
an insulating barrier between the normal and supercon- 
ducting regions. We find that also the differential shot 
noise exhibits a reflectionlcss-tunneling effect observed as 
an enhancement of the noise at zero voltage. Our calcu- 
lations are consistent with the quasiclassical results and 
other previous work. 
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